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Recent theoretical work which employs the classical magneto-ionic theory for a special 
model of the ionosphere applicable to transmission via the ionosphere at or close to grazing 
incidence is employed to analyze LF propagation data. The results of the analysis i 1 lust nil e 
a practical model of the ionosphere by a detailed study of transmission via the first time- 
mode in particular. 

1. Introduction 

Recent contributions and extensions to the classical geometric-optical theory of propaga- 
tion [Wait and Conda, 1958; Wait, 1900a; Johler and Walters, 1960, 1961] of LP waves about 
the earth provide a valuable analysis tool for the interpretation of various experimental data. 
This paper summarizes the mathematical formulas required to predict the field with particular 
emphasis on a type analysis most suitable for large scale electronic computers. Attention is 
focused on certain available experimental data and the results of the analysis are detailed with 
emphasis on technique. However, the physical phenomena which can be deduced as a result 
of the analysis of these data are given considerable attention to emphasize the value of the 
technique. 

2. Mathematical Theory 

A geometric-optical theory [Bremmer, 1949] can be employed with the aid of a sharp]} 
bounded model ionosphere [Johler and Walters, I960, 1961] emplaced at various altitudes h 
above a spherical earth to describe the cw field, E(oo,d), radiated from a dipole source current 
moment, I L It is common engineering practice at LF to relate the amplitude of the field, 

\E(u,d)\ to the well known concept of radiated power [Ballentine, 1924], neglecting earth 
losses 

P r =1.6(10- 13 )c 2 W) 2 /Z , (1) 

where Z ^120ir } a constant. Such a field, E(<a,d), implies a similar Hertzian dipole receiver 
(vertical polarization). The total field, E(a,d), is the sum of j-\-l reflected ionospheric waves, 
j=0, 1, 2, 3 . . ., 

E^D^Ej^d), (2) 

in which the zero order (j=0) reflection, E (u,d) is the groundwave. Each skywave reflection, 
i=l, 2, 3 . . ., using the positive time function, exp (iwt), can be written, 

Ejfad) =i<adDj 1 exp (icd',) G)G] a^FjCj, (3) 

where [Johler, Kellar, and Walters, 1956], 

E (u,d) = [E Q (u,d)} exp [-i(ub+cj> c )], (4) 

and where, 

C=I lb 2 /4wKd*= (10- 7 )/d,(U= 1) • (5) 

The local time x for the groundwave, to is, 

t^t-b, (6) 

b=md/c, (7) 

*This is a revised version of a paper presented at the Avionics Panel Meeting of the Advisory Goup for Aeronautical Research and Develop- 
ment (AGARD) NATO, Istanbul, Turkey, October 1960. 

i The notion of local time was apparently introduced by Lorentz [1906, p. 57]. 
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Figure 1. Relative sky wave delay for various altitude, h, of the lower boundary of the model ionosphere. 

where c is the speed of light, c~3(10 8 ) m/sec, and 771 is the index of refraction of air, iji~l. 
Similarly, the local skywave time, or skywave delay, £,'(j=l, 2, 3 . . .) is, 



tj=t — bj, 
b^nDJc, 0=1,2,3 . . .)• 



(8) 
(9) 



The quantity, bj—b, is frequently called the relative skywave delay, figure 1 (relative to the 
groundwave). The physical length of the ray, D h can be evaluated geometrically for a reflec- 
tion at an altitude, h, above the surface of the earth of radius, a, 

Dj=2j [{a-\-K) cos 4> it j—a cos rj, (10) 

where, figures 2, 3, 4> t is the angle of incidence of the "ray" on the ionosphere and tj is the 




20C 300 400 500 600 TO 

OISTANCE ,d/j , STATUTE MILES 



Figure 2. Geometric-optical relation between the 
angle of incidence, 4> if distance from the source, 
d/j (j = order of sky wave hop or time-mode, ] = 1, 2, 
3 . . . ) and altitude of the boundary of the model 
ionosphere, h. 




Figure 3. Coordinate systems at the ionosphere 
boundary. 
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corresponding angle of incidence on the earth and the subscript "j" reminds the reader that the 
cquation refers to the particular ionospheric reflection under consideration. The angles 
<t>i and tj are evaluated quite simply from the geometry, 

n 

sin ifi =A7 1 a sin — .> (11) 

cos *i li =A7 1 [a(l— cos 0/2j)+A], (12) 

sin T j =Ay 1 (a+h) sin 0/2j, (13) 

cos r^=A7 2 [a(cos 6/2j—l)+h cos 0/2j), (14) 

A,= [2a(a+A)U-cos d/2j)+h 2 ]*, (15) 

where is the angle at the center of the earth subtended by the distance, d, along the surface 
of a spherical earth, or, simply, 

d=ad. 

The factors, G) and G) refer to the transmitter and receiver complex antenna patterns re- 
spectively, and again the subscript, "j", is employed to designate the particular time-mode 
under consideration. The vertically polarized Hertzian dipoles considered in this paper reduce 
to, Gj=G r j=sm Tj. Of course, other antenna complex pattern factors can be introduced. 

Since plane reflection coefficients are employed in this analysis, the convergence-divergence 
coefficient, a jf is necessary to convert the plane reflection coefficients to spherical reflection 
coefficients or from the viewpoint of "ray theory" it is necessary to take account of ray focusing 
by the ionosphere and a corresponding defocusing of rays by the earth. The "classical" con- 
vergence-divergence is modified by a " convergence correction" [Wait, 1959], A h in this paper, 
since the values of the first time-mode studied are at distances close to the "caustic" of rays or 
geometric-optical horizon, and hence the complete expression, a h 

aj= (1 +h/a)[(2j sin 0/2j)/sin 0]*X { [a(l -cos 6/2j) +h]/[(a+h) cos d/2j-a] } *A U (16) 

in which Aj can be evaluated from the cylindrical Hankel function of order, n=\, of the second 
kind, 

Aj^ \ - fact cos 3 tj/3 sin 2 tj I Hi 2) k\d cos 3 r ; /3 sin 2 t } 

Xexp -j -i |5ir/12— fee cos 3 r,/3 sin 2 r/1 V (17) 

The factor, ] A h approaches unity (A^X) at shorter distances, (tj^t/2). Near the geometric- 
optical-horizon, r^7r/2, and beyond, the correction becomes quite appreciable. The Hankel 
functions, H^Hz) can be readily evaluated on an electronic computer by a consideration of the 
integral form, 

If 7 " If 00 
H { n\z)=- exp (— iz sin p+inp)dp z exp (— z sinh p){exp wp + exp (—np+imr)}dp, 

IT Jo Tl j o 

(18) 

where the order, n=$, and the argument, z=k\a cos 3 r^/3 sin 2 r h in which the wave number, 
ki=G)T}i/c^cA)/c, for the case under consideration are real numbers. Nevertheless, it is inter- 
esting to note that the method is also applicable to complex argument, z, and complex order, 
n, provided the multiple branches of the Hankel function are carefully considered. The follow- 
ing substitution is made, 

p=exp (-p), 
and a resultant expression for the Hankel function, 
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H®(z)=± £ exp (-iz sin p+inp) dp-X £ [V»-i+|8-i exp (m*-)] exp |"-| 0-/?)] dft 

(19) 
or, 

#f(s)= f T f(z,P,n)d P - I \f(z,p,n)dfi, (20) 

Jo Jo 

is found to comprise integrals with finite limits. These integrals can then be evaluated in terms 
of Gaussian quadrature, 2 

M M 

H™(z) = T, W^[Uef(z,p n ,n)+i Im f(s, Pm ,n)]-j:Wi[B.ef(z,fi m ,n)+iImJ(z,^,n)]+e(M) (21) 

m=l m=l 

m=l, 2, 3 . . . My where, e(M) can be made arbitrarily small by increasing M, 

where, 

WSr=wHj2 t 

Wl=HJ2, 

p m =wxj2+T/2, 
P m =xj2+l/2, 

and the constants, x m , of the theory of Gaussian quadrature can be evaluated as the roots of 
the Legendre polynomials defined by, 

°~ (x*-l) M -2 M M\P M (x)=0, (22) 

or, 

Po(*)=l, 

3 1 

P 2 ( X )=-X 2 --y 
Ps(x)=^X S —-Xy 



P 4(X)== | X 4_^ 2 + | 



(23) 



where polynomials of higher degree are determined by the use of the recursion formula, 

(M+l)P M+1 (x)+MP M . 1 (x)-(2M+l)xP M (x)=0. (24) 

The weight coefficients, H m , are evaluated from the roots, x m 

H m =2/(l-xl) [P' M (x m )]\ (25) 

The factor, F j9 accounts for the presence of the earth at the transmitter (source) and re- 
ceiver. If it can be assumed that the "ray" reflected by the earth is not too close or beyond 
the geometric-optical horizon, the Fresnel approximation of the ground reflection coefficient 
suffices to determine F Jf or, 

F s ~[1+RI(t)][1+R;(t,)], (26) 



'♦This type analysis has been used extensively by the author in previous papers to evaluate Fourier integrals [Johler and Walters, 1959], 
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where the superscripts, t and r, refer to the transmitter and receiver respectively and the 
subscript, e, refers to vertical polarization, and, 

R e (rj)= [k\ cos Tj/kl-lkl/kl-sm* r^}/{k 2 2 cos r#?+[£|/^--sm 2 77]*} , (27) 

or for horizontal polarization, 

' R m ( Tj )= {cos Tj-lkj/kt-sin 2 r y ]*}/{cos r y +[^/^-sin 2 r,]*} , (28) 

where, 



(29) 



and e 2 is the dielectric constant of the earth (e 2 ~15) and a is the conductivity of the earth 
(<r~0.005 for land). Since the examples considered in this paper involve distances which, 
for the first time-mode, were close to the geometric-optical horizon, it was found necessary to 
replace the Fresnel reflection coefficients with a more rigorous treatment. 

Close to the geometric-optical horizon a calculation of F J =F t j Fj can be accomplished by a 
numerical evaluation of the contour integral, [Wait and Conda, 1958], 



Fy r ~(w)-i exp [-ikiad'] \° exp [-i(ha/2)Wfi]/W{(p)^^Vi(p) dp, 

J 00 exp [ — i'lrr/'d] 



(30) 



d'=(d— d H )ja, where d H is the distance 3 from transmitter to the geometric-optical horizon 
(tj^tt/2) and, 

./kxd\iki I k\ 
Z=-\-2) KV 1 ""** (31) 

in which <x=a 2 is the ground conductivity at either transmitter F) or receiver F] and e 2 is the 
corresponding dielectric constants (relative to a vacuum), and, 

Wi(p)=exp [-2rt/3] V^/3 (-p)Wf [2/3(- P )*], (32) 

W[(p)=exp [-2^/8] VW3{-^(-p)-Wf JJ(- p )f] + p//;« g ( _ p) }J|, 

W[(p)=ex V [-ibc/S] V^73 pi7| 2) [| (- P )»] (33) 

H^(Z)=-^ fff (Z)=exp [-i2x/3]fff (Z)-^ fff (Z). (34) 

The Hankel functions H^ h2) (Z) can be calculated by previously described quadrature methods- 
(20-27). The contour integral can be written [Wait, 1960b], 

mt r / n i r -j nn f f m exp [- i(*ia/2)*0'a] , f OT exp [-t(£,a/2)*0V] , , \ 

(35) 

a'=a exp (— i27r/3). These integrals can then be evaluated with a Gaussian quadrature. 

Close to but beyond the geometric-optical horizon (tj^tt/2) [Wait and Conda, 1958], a 
calculation can be performed as a residue series summation, 

i- ^% exp [— ^(fcicOfo'rf] , v 

i* 7 ;- ' 2^ Vtt exp [-tfc(d-^)] Zgf (2* r f- 2 2)Wi(2*r?) ( 



3 a ~ 6.367 (10 6 ) m. In diffraction problems it is common practice to account for the vertical lapse of the permittivity of air at the earth's sur- 
face by use of an "effective earth radius", a e in place of a such that — ~0. 75— 0.85, [Bremmer, 1949, p. 145]. 
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where the complex conjugate t s of r* comprises the special roots of the differential equation 
of Riccati tabulated by Johler, Walters, and Lilley [1959], employing the time function, 
exp {—i(d), 



where 

8 e = 



-2r s d 2 e +l = 0, (37) 

(IT S 

%k\a^k\ . Z 



*«>'[f-] ! <w'[3£] ! <38) 

and, ki=--7]i^-,k2i implies the complex conjugate of k 2 , (3), and the limiting roots are found 
c c 

from the equations, 

Hf [|(-2r,)*]=0,5 e =c», 

fff[|(~2r s )l]=0,5 e =0, 

s=0, 1, 2, 3 ... . 

Finally, the effective reflection coefficient of geometric-optics, C h can be evaluated for 
vertical polarization, 

U > j\B e dx±l-A 2 x-A^_\ x J w "> 

where, 

a — j? f 

-^J-l — Xl m- 1 - mm 
jri2 == -tt e J. ee\At m l mm 

J±z — ii e £i m l 1 ee l mm -\- A em l me \. (40) 

Thus, the essential nature of the propagation about the earth via the ionosphere can be described 
in terms of four reflection coefficients, T ee , T em , T me , and T mm . The reflection coefficient, T cei 
refers to the vertical electric polarization of the incident plane wave and a similar vertical 
electric polarization of the reflected wave. The coefficient, T em describes the generation of the 
abnormal component by the incident vertical polarization (reflected horizontal polarization 
for vertical excitation). Similarly, T mm refers to the incident horizontal electric polarization 
and the corresponding reflected horizontal electric polarization. Also, the abnormal component 
generated by the horizontal electric polarization (reflected vertical polarization for horizontal 
excitation) is described by the coefficient, T me . Thus, figure 3, 

J- ee—-LLy'r/-fc>y'i } ^m-^'rPi'i, J- em = -fcx'r/&y'i, A me — &y' rl -tLx' i (41) 

where the subscripts i and r refer to incident and reflected wave respectively at the ionosphere 
boundary. The calculation of these reflection coefficients has been previously described in 
detail and the results of various computations have been tabulated [Johler and Walters, I960]. 
The reflection coefficients are completely specified by the angle of incidence <j> iy figures 2, 3, 
electron density iV(El/cm 3 ), figure 4, collision frequency v(c/s), figure 5, magnetic field intensity 
H m (gauss), magnetic dip / and magnetic azimuth <j> a . The effective reflection coefficients, 
Cj, can then be written in terms of the ionosphere reflection coefficients, T, and the ground 
reflection coefficients, 7?, for each local ionosphere or ground reflection, 

Ci=Ci{Ti) = T ee} 

C 2 = 62 (t 2 ) = l eei l ee 2 ^e 1 I ^m x A em Y A me 2 » 

^3 ^3 V 3/ -^^ m l ■*- ee%-*- em± -*■ me^T t\ ^tn2 ee l em l me^l AX- e ^li e ^ A ee ^ 1 ee 2 -*- ee% 1" Alm^ti m ^l m^ emj A me%j 

(42) 
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Figure 4. Electron density-altitude, N(h), profiles 
of the lower ionosphere deduced from various radio 
and rocket measurements. 



NICOLET (1958) 




J L 



J I }sJ_ 



80 90 

ALTITUDE, h, km 



Figure 5. Collision frequency-altitude, v{h), profiles 
of the lower ionosphere deduced from theory and 
various rocket measurements. 
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where the subscripts 1, 2, 3, 4 . . . refer to the first, second, third, fourth, etc., reflection 
regions of the ionosphere with corresponding separate values for electron density, N, collision 
frequency, v, and the intensity and direction of the earth's magnetic field vector, H m or, for 
the ground, refer to the separate values of conductivity, a, and dielectric constant, e 2 . 

For a homogeneous ionosphere and ground with a single set of values, A T , v, H m , <r, and 
e 2 , the effective reflection coefficients, C h reduce as follows: 

C\ = T ee 

1^2 -*- eette\lirn-*- em-*- me 

p 9ft ft T T T -Lft^T^ 4-ft 2 T T T 

^3 ^J-^e^^m-*- ee-*- em-*- rael 1 ^ 1 ee\ - Li, m- L mm-*- em-*- mej 



p ft ft2 T2 nr2 _|_ OD2D nr<2 'T rp A_OftT? 2 T T T T J_P3 7^2 f T 1 I 7?37»4 

^4 ±i/ e Li >m ,*- em J -me\^* J±t e ±i m- L ee-*-em J -me I ^ ±i e ll m J -ee- L mm J -em J -me\ 1 * / m J - mm-*- em ± me\ ±i e J - ee 

ii 

The wave number of the ionosphere propagation medium, 

CO 



(44) 



(45) 



implies two upgoing waves excited in the model ionosphere with index of refraction, jj , for 
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the "ordinary" wave and index of refraction, r] e , for the "extraordinary" wave. The value 
of the reflection coefficients, T, are related to the rj and r\ e . 

The rate of attenuation and phase lag of the wave, E t , transmitted into the ionosphere 
can be deduced by an interpretation of the transmitted wave in terms of the complex index 
of refraction, 77 = 770,,,' 

E t = \E t \ exp / i [ut-2 nDJy. (46) 

The amplitude is as follows: 

exp < - Im rjD >, exp {nepers}, (47) 

or the amplitude change is, 

exp -< - Tm 77 10 3 >> exp {nepers/km}, (48) 

^ c J 

or, 8.686 - Im rjD, decibels, (49) 

or, 8.686 - Im 77 10 3 , db/km, (50) 



where D (meters) is the penetration distance into the ionosphere. The corresponding phase 
lag is, 

CO 

c 
or, 

CO 



Re rjD radians, (51) 



c 



Re 77(10 3 ) radians/km. (52) 



3. Analysis of Experimental Data 

Although there is a tendency in LF propagation studies to separate arctic-auroral phe- 
nomena from those phenomena observed at LF in temperate latitudes, the mechanism of propaga- 
tion seems to differ only in the degree and frequency of the so-called disturbed and blackout 
conditions which characterize the arctic regions. Certain data on propagation at LF in both 
the arctic-auroral regions and the temperate regions are analyzed in this paper to demonstrate 
principle rather than deduce comprehensive engineering conclusions. A theoretical investigation 
and analysis is therefore attempted to gain some insight into LF propagation phenomena for 
Loran-C system evaluation. 

Ionosphere electron density-altitude profiles N(h), figure 4, and collision frequency-altitude 
profiles v(h), figure 5, have recently been deduced with the aid of a theoretical method and an 
experimental HF signal transmitted through the ionosphere from rockets. Waynick [1957] 
reported the "Houston composite/' figure 4 N(h) profile, and Seddon and Jackson [1958] have 
reported "NRL 50" and "NRL 47" profiles. Seddon and Jackson have also reported 
"NN3.08F" v(h) profile, figure 5. Other v(K) data are shown from Crompton et al. [1953], 
Gardner and Pawsey [1953], Fejer [1955], and Nicolet [1958]. 

The "Houston composite" N(h) profile implies an interesting low electron density region 
at the 65 to 70 km altitudes of the lower ionosphere with values between N=10 and N= 100 
El/cm 3 , figure 4. Such a profile seemed to persist in the propagation data examined and pro- 
foundly influenced the daytime ionospheric propagation. However, it does not appear to exist 
in the "NRL 50" data which seems to give rise to great electron densities quite abruptly. 
Unfortunately, as a result of rocket drag considerations, the HF antennas were not extended 
until the rocket attained an altitude of about 55 km, and measurements of electron densities 
in this region are not very detailed. This lack of data below 60 km is more apparent in the 
"NRL 47" data, observed in the auroral zone at Churchill, Canada, which illustrate a very large 
enhancement of the electron density during auroral "blackout" conditions. A "mean of the 
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composite and blackout" profiles was employed as an intermediate "disturbed" model 
ionosphere. 

Various collision frequency-altitude profiles, v(h), are shown in figure 5. The "Nicolet/3" 
curve was used for all model ionospheres employed in this analysis. This is a theoretical curve 
substantiated with experimental rocket data. This theoretical curve is also consistent with 
other findings [Ratcliffe, 1959]. 

Theoretical field intensities, figures 6 to 9, were calculated for comparison with available 
data 4 on LF transmission from Adak to Nome and Kodiak, Alaska. The independent variable 




ADAK-KODIAK 

QUIET 

V= I AMPERE - METER 

WITH FOCUSING CORRECTION 

WITHOUT FOCUSING CORRECTION 



— 







ADAK-KODIAK 
QUIET 

WITH FOCUSING CORRECTION 
• WITHOUT FOCUSING CORRECTION I 




ALTITUDE, h, km 



70 75 80 

ALTITUDE, h, km 



Figure 6. Total theoretical field E(co, d) for various fixed heights, h, of sharply bounded model ionosphere for 

LF transmissions from Adak to Kodiak, Alaska. 

(See figures 4 and 5 for electron density N . . . .) 



« The data employed in this analysis is treated somewhat differently from previous authors. The authors wish to thank A. D. Watt and 
J. E. Bickel for the use of their original data. See also Bickel [1959] and Watt et al. [1959]. 
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of the calculations was the altitude, h, of the assumed sharply bounded model ionosphere. 
This implies a variation of the electron density according to the N(h) profile as the variable 
h is changed. Also, small adjustments in the angle of incidence, <fo, are implied. 

Geophysical data on the ionosphere employed in the calculation are shown, tables 1, 2. 
The magnetic data, shown in table 2, were scaled from Magnetic Charts [1954, 1957]. Appar- 
ently a nominal intensity value H m —Q.b gauss (table 2) is an adequate assumption for many 
practical cases. 

Table 1 



h 


V 


2V, 


N d 


N b 


km 
55 


c/s,X10+6 
175 
50 
24 
16 
10.9 
9.2 
7 

4.5 
2.6 
1.6 
1.2 
0.88 
.465 
.155 


E1/C77&3 


El/cr»3 


El/cm* 

10 


60 






1000 


65 
67.5 


10 

56 

150 

180 

200 

160 

140 

500 

3500 

11000 

15000 

10000 


1205 


2400 


70 
71 


2675 


5200 


72.5 






75 
77.5 


4505 


8850 


80 
81 


7950 
10550 
16100 
20500 


15400 


82.5 




85 
90 


26000 









Table 2. — Magnetic field, H m , data 
Adak-Nome, rf=965.5 statute miles (1,550 km) 





Reflection 
number 


H m 


/ 


<t>a 


1 

2 

2 

9 


1st 

1st 
2d 
1st 
2d 
3d 


Gauss 
0. 5187 
.4996 
.5329 
.4952 
.5187 
.5388 


deg 
68.68 
65.75 
71.52 
64.82 
68.68 
72.46 


deg 
12.27 
12.52 
11.72 
12.63 


3... 

3 


12.27 

11.66 



Adak-Kodiak, <2 = 1040.1 statute miles (1,670 km) 



1 

2 

9 


1st 
1st 
2d 
1st 
2d 
3d 


0. 5035 
.4950 
.5170 
.4923 
.5035 
.5262 


67.18 
65.15 
69.04 
64.42 
67.18 
69.64 


51.08 
52.84 
49.73 


3 

3 


53. 42 
51.08 
49.35 







The modal (skywave time-mode) and the total (including the groundwave) calculated 
fields (amplitude and phase) for quiescent, disturbed, and blackout conditions of the ionosphere 
are presented, figures 6 to 9. The effect of the ionosphere focusing corrections, A h is also 
illustrated, and, indeed, was found to be appreciable, both for the total field, E(ca,d), and the 
first modal field, Ej(u,d),j = l. The diffraction correction i^ was also modified from the clas- 
sical Fresnel reflection coefficient as discussed with the aid of calculated data of Wait and Conda 
[1958] which was installed in the electronic memory of a small computer. 5 

The ground wave, E (u,d), which was calculated with the classical Bremmer-van der 
Pol theory [Johler, Kellar, and Walters, 1956] and corrected for land-sea boundaries with the 
aid of the Millington [1949] method is not shown alone graphically, since the groundwave is 
not a function of the emplacement altitude, h r but can be summarized as follows: 

Table 3 





E n (u>,d) 


arg E (o>,d) 




Iol = l 


7oZ = 2050 




Nome 


v/m 

7.56(10-9) 
4.64(10-9) 


v/m 
1 .55(10-5) 
9.51(10-6) 


Radians 

2.30 


Kodiak. 


2.43 



5 IBM-650. 
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All computations were made with the assumption of a dipole current moment, I l=l. 
The total radiated power, Pr =1.3 kw, for the arctic data considered was found to correspond 
to a dipole current moment, I l=2050 amp-m; therefore, the calculated field amplitude, [E(co,d)], 
or [Ej((a,d)]j has been increased by this factor for comparison with the data. 
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Figure 7. Disturbed-blackout theoretical field E(w, d) for various fixed heights, h, of sharply bounded model 
ionosphere for LF transmissions from Adak to Kodiak, Alaska. 



517 



E 

^ -9 

> icr 




AOAK - NOME 

QUIET 

I / : ' AMPERE -METER 

WITH FOCUSING CORRECTION 

WITHOUT FOCUSING CORRECTION 



65 70 75 80 85 90 

ALTITUDE, h, km 



A DAK -NOME 

QUIET 
-— WITH FOCUSING CORRECTION 
WITHOUT FOCUSING CORRECTION 




65 70 75 80. 85 90 

ALTITUDE, h, km 



Figure 8. Total theoretical field E(co, d) for various fixed heights, h, of sharply bounded model ionosphere for LF 

transmissions from Adak to Kodiak, Alaska. 
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Figure 9. Disturbed-blackout theoretical field E(<o, d) for various fixed heights, h, of sharply bounded model 
ionosphere for LF transmissions from Adah to Kodiak, Alaska. 
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Figure 10. Comparison of quiescent and disturbed-blackout conditions of theoretical reflection coefficient, T, 
(j = l), for various fixed heights, h, of sharply bounded model ionosphere for LF transmissions from Adak to 
Kodiak, Alaska. 

H m =0.5035, a =51.O8°, 7=67.18° 
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In addition to the total field, E(co,d), and the modal field, Ej(co,d), figures 6 to 9, the 
reflection coefficients, T, (horizontal and vertical polarization; normal and abnormal com- 
ponents), figures 10 to 11, and the transmission data of the ordinary and extraordinary upgoing 
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Figure 11. Comparison of quiescent and disturbed-blackout conditions of theoretical reflection coefficient, T, 
(] = \), for various fixed heights, h, of sharply bounded model ionosphere for LF transmissions from Adak to 
Nome, Alaska. 
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Figure 12. Attenuation rate and phase lag rate at various emplacement altitudes, h, for upgoing "ordinary" and 
"extraordinary" waves excited in the model ionosphere by Adak to Kodiak, Alaska transmissions during 
quiescent and disturbed-blackout conditions. 
First time-mode: H m =0.5035, 0«=51.O8°, 7=67.18°. 
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Figure 13. Attenuation rate and phase lag rate at various emplacement altitudes, h, for upgoing "ordinary" and 
(l extraordinary 7 , waves excited in the model ionosphere by Adak to Nome, Alaska transmissions during quiescent 
and disturbed-blackout conditions. 
First time-mode: H m =0.5187, O =12.27°, 7=68.68°. 
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waves excited in the ionosphere by the reflection process, figures 12, 13, were evaluated. The 
unabsorbed modal field, Ej(oo,d)/C j} and the effective reflection coefficients, C h are separately 
illustrated, figures 14 and 15, respectively. 

The emplacement height, h, of the model ionosphere was deduced by comparison of the 
computed field strength, E(u,d), with the measured value. The corresponding four reflection 
coefficients, T ee , T em , T mm , T me , (amplitude and phase) were then evaluated from the theory. 
The results of this process which were readily scaled from graphs such as illustrated, figures 10, 
11, are presented in table 4, and the reflection coefficients deduced for the first time-mode only 
are presented. 

The higher order time-mode, .7—2, 3, reflection coefficients are obviously quite numerous 
since a separate set exists for each ionosphere reflecting region. These reflection coefficients 
are determined by the analysis by reference to the emplacement altitude, h, table 5, deduced 
from the field, E(u),d). 
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Figure 14. Unabsorbed modal field, Ej(o>, d)/c, at 
various model ionosphere altitude emplacements, h, 
for Adah to Nome and Kodiak, Alaska transmissions 
\l=l ampere-meter. 



Table 4 



Adak-Kodiak (Daytime) 



%Time 
Field \E\ 
exceeded 


Season and year 


Observed 
|£(«,d)| 


h 


arg E{u!,d) 


1 #i(o>, d) | 


arg #i(«, d) 


\Tee\ 


arg Tee 


\T em \ 


arg Tem 


\T mm \ 


arg T mm 


\Tme\ 


arg T mt 


10 

50 

90 

10 


Winter 1953-54.. 
Autumn 1954 


fiVlm 
33 
12 

4.8 

20 
7.1 
2.6 


km 

69 

68 

<65 

69 
<65 
<65 


Radians 
10.4 
9.5 
<5.0 

10.4 
<5.0 
<5.0 


nvlm 
56 
35 
<2. 

56 

<2.0 
<2.0 


Radians 
2.8 
2.6 
<2.1 

2.8 
<2.1 
<2.1 


0.27 

.19 

<.02 

.27 
<.02 
<.02 


2.1 

2.0 

<1.6 

2.1 

<1. 6 
<1.6 


0.03 
.02 


.03 




5.8 
5.7 

<4.8 

5.8 
<4.8 
<4.8 


0.23 

.16 
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.23 
<.02 

<.02 
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<1.6 
<1.6 


0. 03 
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<1. 6 







Adak-Nome (Daytime) 



10.. 
50-. 
90.. 

10.. 
50.. 
90_. 

10_. 
50-. 
90_. 



Winter 1953-54.. 


124 


69 


9.8 


94 


3.0 


0.29 


2.1 


0.03 


5.95 


0.23 


1.95 


0.055 




56 


68 


9.0 


60 


2.7 


.18 


1.9 


.02 


5.55 


.15 


1.8 


.04 




22 


65.5 


6.4 


9.0 


2.4 


.04 


1.6 


.01 


4.8 


.035 


1.6 


<.01 


Autumn 1954 


64 


68. 5 


9.3 


70 


2.8 


.22 


1.95 


.025 


5.7 


.18 


1.9 


.045 




28 


66.5 


7.4 


22 


2.5 


.08 


1.65 


.01 


5.0 


.065 


1.65 


.015 




11 


<65 


<5. 7 


<5.0 


<2.3 


<.02 


<1. 55 





<4. 75 


<.02 


<1. 6 





Spring 1954 


34 


67 


7.8 


36 


2.6 


. 11 


1.75 


.015 


5.2 


.09 


1.7 


.02 




19 


<65 


<5.7 


<5.0 


<2.3 


<.02 


<1. 55 





<4. 75 


<.02 


<1.6 







10 


<65 


<5.7 


<5.0 


<2.3 


<.02 


<1. 55 





<4. 75 


<.02 


<1.6 






2.25 
2.0 
1.6 

2.1 
1.75 
<1.6 

1.8 
<1.6 
<1.6 



523 



1 r 



i i i r 




ADAK-KODIAK 
1 



180 


1 


1 1 1 


1 1 


16.0 


- 




- 


14.0 






f ^^ j = 3 ~~ 


12.0 


- 




- 


100 






- 


8.0 






X* *"** j=2 _ 


6.0 






- 






^r 


^ j«l ~~ 


40 




2.0 

no 


1 


1 1 1 


ADAK-KODIAK 

1 1 



65 70 75 

ALTITUDE, h, km 



65 70 75 80 

ALTITUDE, h, km 




5 12 



i r 




J I I I I L 



65 70 75 80 85 

ALTITUDE, h, km 



Figure 15. Effective reflection coefficient, Cj, for Adak to Nome and Kodiak, Alaska, transmissions under 

quiescent ionosphere conditions. 
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Table 5 



Condition 


Mode 


N 


h 


|E(«,d)| 


arg E 


b 


Total _._ 
Total .._ 

3=1 

Total. ._ 

i=i 


El/cm3 
10 
1000 
10 
10 
10 


km 
55 
60 
55 
65 
65 


v/m 
7. 59 (10-9) 
6. 60 (10-*) 
1. 18 (lO-io) 
9. 77 (10-9) 
2. 42 (10-9) 


radians 
5.54 


b 

b 

q__ .__ _ 


2.84 
2.10 
5.66 


q ... _. _ 


2.33 











It is interesting to note emplacement altitudes, h, less than 65 km under quiescent conditions. 
Indeed, the highest altitude for these daytime values was 69 km (winter 1954). A good nominal 
daytime emplacement altitude is 67.5 km. The low value (|T|<0.1) reflection coefficients are 
also quite interesting. Note the highest value, \T\~0.27. Note also, the median value at 
Nome (1954 equinox), |T ec | = 0.08. A low value |7 7 |<0.1, and "black", \T\~0 appearance of the 
ionosphere indicates a low level absorbing region in the lower ionosphere at grazing incidence. 
A region of low electron density, N= 10-100 El/cm 3 , is indicated by the model emploj^ed, 
figure 4. The transmission into such a region at grazing incidence, 0^80°, figures 11, 12, 
indicates almost complete absorption of both the "ordinary" and the "extraordinary" rays if a 
"virtual height," h-{-Ah, is assumed of the order of 70 km. The transmission curves, figures 
11, 12, demonstrate that the "extraordinary" ray is very highly attenuated such that it is 
necessary to consider only the "ordinary" ray. Assuming 1 db/km, for example, figures 11, 12, 
and a virtual height, h+Ah, of 70 km, noting that appreciable attenuation rate begins at 65 km, 
the total attenuation can be estimated; total decibels ~2A&/cos </>^57 db. This, of course, 
neglects a reflection absorption at the 70 km virtual height plus the transmission absorption 
back into the region below the ionosphere. The values, Ah, can be decreased and a correspond- 
ing decrease in the attenuation noted. But a repetition of the process leads to the interesting 
conclusion that most of the ray reflection must occur in a fraction of a kilometer, and, hence, 
the sharply bounded model can be reasonably applied under these conditions to predict the LF 
field, E(o),d). A similar examination of the second and third time-modes leads to an almost 
similar reflecting region since the angle of incidence, <j> u figure 2, does not change rapidly until 
the distance, d/j, is sufficiently short. It is, nevertheless, conceivable that the amplitude 4 of the 
skywave reflections, ,7 = 4, 5, 6, . . . at shorter distances could become quite large as a result 
of reflection from a high level "virtual height," h-\-Ah=70 to 85 km. Under these conditions, 
however, the sharply bounded model would not be valid and it would he necessary to develop 
theoretical techniques for a diffusely bounded model ionosphere. 

Information on the lower ionosphere N(h) profile during nighttime conditions is quite 
sparse and hence the theoretician can only speculate on its nature. However, under certain 
reasonable assumptions, a possible explanation of the increase in field strength experienced at 
night may be deduced from the daytime model ionosphere employed. 

Despite the seeming "evaporation" of the lower regions of the ionosphere at night with a 
corresponding abrupt increase of the electron density to a large finite value, it seems reasonable 
to assume that this abrupt increase does not occur within A/10 or 0.22 km at 135.6 kc/s. This 
makes it possible to utilize the present sharply bounded model ionosphere, if il is emplaced at a 
higher level. Recent work [Ellyett and Watts, 1959] indicates a height of 80 to 85 km is not 
unreasonable. 

Since the change of the angle of incidence, 4>u figure 2, with ionosphere height, h, is not 
great at this range, we may neglect its effect on the field strength, figures 6 to 9, and the effective 
reflection coefficients, figure 15. The only factor changed, therefore, is the un absorbed field, 
figure 14. It therefore seems reasonable to expect the first time mode, \Ei(<a,d)\, to increase 
on the order of at least twice, whereas, reference to the arctic data and figure 16, "Battle Lake 
data" (to be discussed below), indicates an increase of perhaps 4 to 10 times. This could easily 
be explained by a slight change of slope in the N(h) profile employed as evidenced by the steep 
slope of the d curves = 1), figure 15, in the region involved. On this basis, it is reasonable 
to expect a nighttime reflection coefficient, ITI^O.3, 0.4. Indeed, at temperate latitude, these 
values are confirmed [Ross, 1959]. 
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It is not surprising that a large phase change could be expected during the sunset-sunrise 
transition periods, figures 6 and 8. Indeed, such phase changes were described at 100 kc/s 
by Doherty [1957], while operating the Loran-C on the first time-mode of the sky wave. 

The behavior of the LF propagation during disturbed and blackout conditions which are 
especially frequent in the arctic-auroral region is not quite so clear. Figures 7 and 9 illustrate 
disturbed-blackout conditions based on the geophysical measurements of electron density 
profiles, N(h), of the lower ionosphere, figure 4. The corresponding reflection coefficients, 
figures 10, 11, and the transmission of upgoing ordinary and extraordinary waves excited by 
the reflection process are also illustrated, figures 12, 13. The disturbance of the earth's 
magnetic field during extreme blackout conditions is negligible since the total field, H m ^0.5 
gauss and the fluctuations are of the order of 0.001 gauss [Chapman, 1940]. The most 
significant change in the propagation is caused by the drastic change in the electron density 
profile, N(h), figure 4. Unfortunately, Mie shape of the lower levels, <60 km, of the "black- 
out" N(h) profile is unknown. Since the analysis under quiescent conditions indicates low- 
electron density reflections, the curve, figure 4, was extrapolated, and the corresponding field, 
E(a),d), calculated for Adak to Nome transmission. The results are tabulated in table 5. 

Note the j=l mode for a region of iV=10 El/cm 3 has undergone a 26 db decrease and 
the total field is therefore almost entirely groundwave. Appreciable as the decrease is, it is 
still well within the range of observations of arctic data. Notice, figure 16 "Battle Lake 
data" (to be discussed below), that the range of observations in temperate latitudes is approx- 
imately 18 db. The additional fluctuation may then, reasonably, be attributed to the change 
in electron density of the extreme lower ionosphere. If this change were to be so abrupt as 
to reach A^=1000 El /cm 3 in, say, X/10 km, the field strength would show a considerable increase, 
instead of decrease. It is reasonable to expect that a more realistic value of somewhat less 
than iV=1000 El/cm 2 would also produce an increase in the received field strength. The 
shape of the N(h) profile, thus, becomes quite important. 
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Figure 16. Amplitude of theoretical total field 
|E(co, d)| as a function of distance, d, from the 
source illustrating a comparison of experimental 
pulse amplitude measurements with theoretical first 
time-mode in the diffraction zone beyond the geo- 
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Figure 17. Phase of theoretical field arg E(co, d) as a 
function of distance. 
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The calculation of the amplitude and phase of the total field, E(<a,d) as a function of 
distance is illustrated in figures 16 and 17. The corresponding reflection coefficient T ee 
(vertical polarization, normal component) is illustrated in figure 18. 

Amplitude measurements 6 during the day which illustrate the diffraction of the first 
time-mode (j=l) near and beyond the geometric-optical horizon are illustrated in figure 16. 
The agreement with the theoretical curve j=l appears to be quite close except at the range 
of 1,800 miles. It is interesting to note the apparent decrease in the range of amplitude 
observations in the daytime, indicating greater amplitude stability with range in the 
diffraction zone. 



heel ■* 
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Figure 18. Amplitude and phase of the normal com- 
ponent reflection coefficient T ee (vertical polarization) 
as a function of distance and frequency for east-to- 
west <}> & = 270 o and west-to-east <j> & — 90° propagation. 
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4. Conclusion 

The analysis of these data demonstrates the application of propagation theory together 
with a suitable model for the ionosphere to the prediction of LF phenomena. The precision 
of the model for the case of transmission via the ionosphere at or close to grazing incidence is 
considered to be adequate upon considering the state of knowledge of the many parameters 
which influence such propagation. Nevertheless, the absence of an abundance of specific 
information does not deter the author from certain general conclusions which can be made as 
a result of this research. 

At oblique incidence on the ionosphere, a condition which is certainly met under most 
practical situations, the reflection of LF waves during the daytime occurs in a region of very 
low electron density and high absorption, with corresponding low values of reflection coefficients, 
ITKO.l. The dominant skywave time-mode is, in general, the first; however, the second, 
third and higher order time-modes may offer appreciable contribution to the total field as the 
angle of incidence decreases as a result of reflection from higher levels and shorter transit 
distance in the ionized medium. 

Periods of ionospheric disturbance (blackout) can produce a decrease in the field as a result 
of a reduction in the amplitude of the ionospheric time-modes. However, it is equally possible 
for the skywave mode to be incremented. This phenomenon is dependent upon the slope of 



3 From the author's persona! notes on the 1953-1954 testing of the Loran-C (Cytac) System. 
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the electron density-altitude profile, N(h) h. However, if ordinary communication is an 
engineering objective, LF seem to be quite reliable since the characteristic high absorption of 
HF waves does not appear to exist. 

Since the writing of this paper the extensive theoretical work of Barron [1961] and the 
extensive experimental and engineering work of Belrose et al., [1959] have come to the attention 
of the author. It is interesting to note quite good agreement, figure 18, between the reflection 
coefficients deduced from the sharply bounded model presented in this paper and those pre- 
sented as a result of Barron's full wave computations. Barron's data does, however, indicate, 
as might be expected, higher attenuation at the higher frequencies (>150 kc/s). 

It should be noted in conclusion that the adequacy of the sharply bounded model is subject 
to severe limitations, especially at the higher frequencies and the results of this paper suggest 
also the introduction into this geometric optical analysis of more rigorous or flexible models 
for the reflection coefficient which would treat a more diffuse ionosphere lower boundary. 
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Selected Abstracts 

Mean electron density variations of the quiet ionosphere 3 — 
May 1959, J. W. Wright, L. R. Wescott, and D. J. Brown, 
NBS Tech. Note 40-8 \PB 15 1899-3) {I960) $1.50. 
The CRPL has initiated a program for large-scale computa- 
tion of electron density profiles from ionospheric vertical 
soundings. Scaling is performed at field stations permitting 
computation of hourly profiles at the Central Laboratory. 
These profiles are combined to form hourly mean quiet pro- 
files for each station and month. The results of this program 
for the month of May are illustrated graphically. This report 
is the third of a series illustrating the electron density varia- 
tions in the mean quiet ionosphere between latitudes 15 °N 
and 50 °N along the 75 °W meridian. 

On the nature of equatorial spread -F,R. Cohen and K. L. 
Bowles, ./. Geophys. Research 66, No. 4, 1081-1106 {April 
1961). 

Ionospheric propagation via scattering from the F region 
was sought at 50 Mc/s over a transequatorial path (with 
midpoint near Huancayo, Peru) employing a 2580-km trans- 
mitter-receiver separation. Propagation via F scatter was 
present over this path about 10 per cent of the time, though 
only at night. A condition closely related to the occurrence 
of F scatter was the presence of equatorial spread-^ configura- 
tions on the Huancayo ionograms. The height of the propa- 
gation mediumsupportingthisFscatter was usually identifiable 
with the lowest height of the associated equatorial spread F 
on the ionograms. On the basis of an interpretation of ex- 
perimental results, equatorial spread F is shown to arise from 
scattering by relatively thin sheets of irregularities in electron 
density which occur at the bottom of the F layer or as much 
as 100 km lower. The thickness of these patches is estimated 
at the order of 50 km. The scattering irregularities comprising 
these patches are found to be elongated along the earth's 
magnetic lines of force. Electron irregularities of scale 10 
meters or smaller measured in at least one dimension trans- 
verse to the magnetic field lines, and 1,000 meters or longer 
measured along the magnetic field lines, are shown to exist 
as high as 450 km in the nocturnal equatorial F region during 
equatorial spread-/' 7 conditions. Equatorial spread-/' 7 echoes 
observed at Huancayo are demonstrated to be arriving at all 
elevation angles in the magnetic east- west plane. The geo- 
graphical extension of a given scattering sheet in the magnetic 
east- west direction was at times as great as 1000 km. An 
observational procedure is suggested for distinguishing two 
fundamental varieties of spread- .F echoes appearing on the 
equatorial ionograms. A necessary condition, that the con- 
tours of mean electron density be parallel to the magnetic 
lines of force, is proposed as a controlling factor for the oc- 
currence of spread F in the equatorial ionosphere. 

On the analysis of polarization rotation recordings of satellite 
radio signals, K. 8. Lawrence and C. G. Little, pp. 391-399 
of Some ionospheric results obtained during the International 
Geophysical Year; Proc. symp. organized by the URSI/AGI 

committee, Brussels, 1959, edited by W. J . G. Beynon {Elsevier 
Publ. Co., Amsterdam, 1960). 

At frequencies as low as 20 mc/s the rotation of the plane of 
polarization of satellite signals cannot be assumed to be pro- 
portional to f N B cos dl along the line of sight. A method 
of analysis is presented which eliminates this assumption 
and permits accurate estimates of total electron content below 
the satellite height to be obtained from 20 mc observations. 
Some of the results obtained using the method to determine 
the sub-satellite electron content, throughout a satellite pass, 
are presented. These analyses permit the study of ionospheric 
tilts, and have revealed the presence of large-scale irregulari- 
ties in the ionospheric electron content. The deviation in 
sub-satellite content of these irregularities is of the order 



2% of the running mean value; the lateral extent of the irregu- 
larities is of the order of a few hundred kilometers. The ef- 
fects of such irregularities upon ionospheric refraction and 
satellite Doppler curves are briefly discussed. 

Prolonged space-wave fadeouts in tropospheric propagation, 

A. P. Barsis and M. E. Johnson, NBS Tech. Note 88 
{PB161589) {1961) $2.00. 

This paper contains the results of studies performed during 
the last several years on the short-term variability of tropo- 
spheric signals received over within-the-horizon paths. Signal 
variations of this type have been termed "prolonged space- 
wave fadeouts," as they are mainly characterized by reduc- 
tions in signal level to many decibels below presumably 
constant values determined from geometrical optics methods. 
The data described here were obtained from measurements 
over propagal ion paths in the Pacific Coast region of Southern 
California, and the continental region of Eastern Colorado. 
Fadeouts are analyzed as a function of carrier frequency, 
path characteristics, and meteorological parameters. The 
study also includes an evaluation of fadeouts observed oyer 
a path using a mountain peak as a diffracting knife-edge like 
obstacle between transmitter and receiver. 

The relationship of low-height ionosonde echoes to auroral- 
zone absorption and VHF D scatter, J. Iv. Olesen and J. W. 
Wright, ./. Geophys. Research 66, No. 4, 1127-1184 {April 
1961). 

Weak diffuse \\F reflections at heights between 75 and 95 
km in the auroral zone are designated E,d. Diurnally, E,d 
occurs most frequently at noon at all seasons, and shows a 
pronounced minimum of occurrence around 2000 hours local 
time. There is little evidence for a seasonal variation from 
the year's data presented here. It is suggested that E s d is 
related to the anomalies of auroral-zone absorption and that 
it is the layer responsible for VHF forward scatter. 

The Ephi system for VLF direction finding, (J. lb thy, 
It. F. Linfield, and T. L. Davis, ./. Research NBS 65C (Eng. 

and Instr.), No. 1, 48~49 {Jan.-Mar. 1961). 
A new system of VLF direction finding has been developed 
and tested by the National Bureau of Standards, Boulder, 
Colorado. The system has been named "Ephi" (E-<£) 
because the bearing of the transient signal is determined from 
the relative phase (</>) of the vertical electrical field (E) 
received at spaced antennas. The advantage of this scheme 
compared to conventional crossed-loop techniques is that it 
minimizes siting and polarization errors. A minimum of 
three antennas must be used to resolve directional ambiguity, 
and the preferred antenna location is at the vertices of an 
equilateral triangle with baselines equal to % to }{q wavelength 
(at 10 kc/s). Appropriate phase detectors, delay lines and 
coincidence circuits are used to obtain a directional code in 
preset sectors. Within practical instrumentation limits any 
number of sectors of variable widths can be operated simul- 
taneously and each can be rotated in azimuth. 

Integrated starlight over the sky, F. E. Roach and L. R. 
Megill, Astrophysical J. 133, No. 1, 228-242 {Jan. 1961). 
Calculations are made of the total integrated starlight over 
the entire sky based on the star counts in Groningen Pub., 
No. 43. The' results are given in both the photographic and 
the visual magnitude scales in tabular and graphical form. 

Resonators for millimeter and submillimeter wavelengths, 

W. Culshaw, IRE Trans. Microwave Theory and Tech. MTT-9, 
No. 2, 135-144 {Mar. 1961). 

Further considerations on the mm-wave Fabry-Perot inter- 
ferometer are presented. Computed Q values for parallel 
metal plate resonators indicate that at spacings around 2.5 
cm, values ranging from 60,000 at 3 mm, to 300,000 at 0.1 mm 
wavelengths are possible. The plates must, however, be 
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quite flat. These results are important for many investiga- 
tions, and in particular for mm and sub-mm wave maser 
research. For the aperture per wavelength ratios possible 
here, diffraction effects should be small. Consideration is 
given to using curved reflectors or focused radiation in appli- 
cations where the fields must be concentrated. For this 
purpose, re-entrant conical spherical resonators are treated 
in detail, as regards operation in the TEM mode at high 
orders of interference. Expressions for the Q and shunt 
impedance are given, and high values are possible at mm and 
sub-mm wavelengths. Quasi-optical methods of coupling 
into and out of such a resonator are proposed, and the higher 
modes possible in such a resonator are considered. Results 
indicate that it could have application to the mm-wave 
generation problem, and that it represents a good resonant 
cavity for solid-state research at mm and sub-mm wave- 
lengths, and for maser applications in particular. 

The electromagnetic fields of a dipole in the presence of a 
thin plasma sheet, J. R. Wait, Appl. Sci. Research, Section B, 
8, 397-4-17 {I960). 

The problem of electric and magnetic dipoles located near a 
thin planar slab or sheet of ionized material is considered. 
A constant and uniform magnetic field is impressed on the 
slab. Under the assumption that the thickness of the slab 
is very small, expressions for the resultant fields are obtained. 
As a result of the anisotropy of the sheet it is indicated that 
the fields are elliptically polarized in general. On carrying 
out a saddle-point evaluation of the integrals in the formal 
solution it is shown that the far fields may be split into 
"radiation" and "surface wave" components. The depend- 
ence of the radiation pattern and the surface wave character- 
istics on electron density, collision frequency and the impressed 
magnetic field is illustrated. 

The Nation's electronic standards program: Where do we 
now stand? Harvev W. Lance, IRF Trans. Instrumentation 
1-9, No. 2, 94-100 (Sept. 1960). 

Recently increased standards activities throughout the nation 
have resulted in a large increase in the calibration services 
requested from the National Bureau of Standards and in a 
number of related problems. Several of these are discussed, 
including desirable criteria for the staffing and instrumenta- 
tion of electronic standards laboratories, the choice of inter- 
laboratory standards to be submitted to NBS, and the use 
of these standards before and after calibration. An attempt 
is made to clarify the concept of "traceability of calibra- 
tions to NBS," and conditions under which a calibration 
may be called traceable are discussed. Suggestions are 
made regarding the use of available calibration services and 
the action to be taken when needed calibration services are 
not available. The need for an association of standards 
laboratories is pointed out, and many areas are noted in 
which such an association could make valuable contributions. 

Propagation of error in a chain of standards, A. G. McNish 
and J. M. Cameron, IRE Trans. Instrumentation 1-9, No. 2, 
101-104 (Sept. 1960). 

There is a widespread belief that in a chain of standards 
traceable back to a central laboratory, each standard must be 
calibrated in terms of a standard 10 times as accurate. 
Conversely, this implies a 10-fold degradation in accuracy in 
each successive calibration. If this were true, accuracy 
would be intolerably degraded at the working level. A more 
realistic view is that accuracy in a chain of standards is as 
good as its weakest link. This can be demonstrated by 
experimental data on the errors obtained in calibrations and 
can be explained by examining the factors which affect 
accuracy in calibration. 

An analysis of the accumulated error in a hierarchy of cali- 
brations, E. L. Crow, IRE Trans. Instrumentation 1-9, No. 2, 
105-114 (Sept. 1960). 

Calibrations of many types are performed in a hierarchy of 
calibration laboratories fanning out from a national standard. 
Often the statement is made that the accuracy of each echelon 
of the hierarchy should be 10 times the accuracy of the immedi- 
ately following echelon. The validity of such statements is 
examined by deriving formulas for the total error accumulated 



over the entire sequence when systematic and random errors 
may occur in each echelon, and by determining how a given 
total error may be achieved at minimum total cost under 
reasonable assumptions for the form of the cost-error functions. 

Atomic beam frequency standards, R. C. Mockler, R. E. 
Beehler, and C. S. Snider, I RE Trans. Instrumentation 1-9, 
No. 2, 120-132 (Sept. 1960). 

A general qualitative description of atomic beam devices is 
given, with particular emphasis upon cesium atomic frequency 
standards. 

The various uncertainties in frequency measurements are dis- 
cussed including the effects of the uniform C field, the effect of 
a phase difference between the two oscillating fields that 
excite the atomic transition, and the effect of exciting an 
atomic transition with signals that are not monochromatic. 
The two cesium beam devices constructed at the National 
Bureau of Standards are described in some detail. The 
various tests are described for determining the uncertainties 
in the frequency measurements. The results of these tests 
are given. 

The accuracy of these rather dissimilar machines is considered 
to be ±1.5X10~ n , and the precision of measurement is 
2X 10~ 12 . The latter figure represents the standard deviation 
of the mean for measurement times of a few hours, or, about 
100 separate measurements. 

The power spectrum and its importance in precise frequency 
measurements, J. A. Barnes and R. C. Mockler, IRE Trans. 
Instrumentation 1-9, No. 2, 149-155 (Sept. 1960). 
The power spectral density functions of a frequency multiplier 
chain, driven by several different crystal oscillators, were 
obtained by comparing the output with a second chain which 
was stabilized with an ammonia maser. The frequency of 
the maser stabilized chain was demonstrated to be relatively 
fixed; the power spectrum of the other chain was determined 
by two different methods. The results are compared. 
Possible errors and uncertainties introduced by the methods 
are discussed. An analysis is made that relates the instanta- 
neous frequency fluctuations of a signal with the power spectral 
density function. 

Analysis predicts that when frequency modulation occurs in 
the first stages of frequency multiplication or in the primary 
frequency oscillator, the output power spectrum is, in general, 
not symmetrical. Furthermore, the sidebands are increased 
in intensity by the multiplication process. This is, in fact, 
observed to be the case. It is shown that a frequency counter 
will measure the frequency of the center of gravity of the 
power spectrum. 

If signals having a complex power spectrum are used in precise 
frequency measurements, errors may result. 

Standards and measurements of microwave surface imped- 
ance, skin depth, conductivity and Q, H. E. Bussey, IRE 
Trans. Instrumentation 1-9, No. 2, 171-175 (Sept. 1960). 
The metal walls of a cavity resonator become standards of 
skin depth when the internal Q is determined, if the walls are 
uniform and if contact loss is eliminated as in the circular 
TEon mode. An end plate standardized in this way may be 
replaced by another flat piece of material and its unknown skin 
depth easily determined. A copper cavity in use has 94 
per cent of theoretical Q. Deoxidization of the surface raises 
the value to 97 per cent. Surface roughness and mode con- 
version may account for the remaining Q departure. Accuracy 
of Q measurement (±3 per cent) will be discussed and pre- 
liminary results on several metals will be given. 

A radio-frequency permittimeter, R. C. Powell and A. L. 
Rasmussen, IRE Trans. Instrumentation 1-9, No. 2, 179-184 
(Sept. 1960). 

A coaxial RF impedance transformer in which the secondary is 
a single turn of the material to be measured is used with two- 
terminal impedance bridges to determine the complex per- 
mittivity or complex conductivity of low-impedance materials. 
No electrodes are needed, and many conductors, semiconduc- 
tors, electrolytes and high-permittivity materials can be 
evaluated to about 1 per cent, since errors due to electrode 
impedance and interaction as well as first-order series induct- 
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ance arc eliminated. The design, calibration, range, and 
accuracy are given along with measured values of such mate- 
rials as ferrites and strong electrolytic solutions, showing 
complex conductivities differing considerably from those 
previously observed by other methods. 

A transfer instrument for the intercomparison of microwave 
power meters, G. F. Engen, IRE Trans. Instrumentation 
1-9, 202-208 (Sept. 1960). 

In the intercomparison of two microwave power standards, 
or in the calibration of a microwave power meter by means 
of a second or ' 'standard" power meter, the measurement of 
microwave impedance has played a major role. 
Through an extension of the reflectometer concept, it is now 
possible to devise a four-arm junction which, when properly 
adjusted, makes possible the intercomparison or calibration 
of such power meters with little or no regard for their im- 
pedance characteristics. In addition, the method is substan- 
tially independent of the impedance discontinuity which may 
be present at the input flange or connector. This latter result 
is of particular value in coaxial systems. 

Absolute measurement of temperatures of microwave noise 
sources, A. J. Estin, C. L. Trembath, J. S. Wells, and W. C. 
Davwitt, IRE Trans. Instrumentation 1-9, 209-218 (Sept. 
1960). 

In recent years, the requirement for noise-figure measure- 
ments of low-noise devices has led to a need for more precise 
knowledge of temperatures of noise sources. A system has 
been developed which is capable of measuring the absolute 
noise temperature of a gas discharge tube in WR-90 wave- 
guide to within 0.01 db. The comparison circuit is a modified 
Dicke radiometer. The standard noise source is a black 
body radiator at approximately 1000 °C. Preliminary meas- 
urements of noise temperatures of a set of nine laboratory 
equivalents of a commonly used argon discharge tube indicate 
an excess noise temperature of 15.90 db at 200 ma discharge 
current with a variation between tubes of 0.03 db. 

A modulated subcarrier technique of measuring microwave 
phase shifts, G. E. Schafer, IRE Trans. Instrumentation 1—9, 
217-219 (Sept. 1960). 

This paper describes a technique for comparing phases of two 
microwave signals which employs amplitude modulation in 
one channel of a two-channel system. This modulated sub- 
carrier technique produces a null response in a two-channel 
system for all ratios of the amplitudes of the waves traversing 
the separate channels. Therefore, high precision is main- 
tained for all ratios of amplitudes. A practical application 
of this technique using readily available components achieved 
a precision of phase measurements of a tenth of a degree for 
a change of attenuation of 50 db of component under test. 

Measurement of reflections and losses of waveguide joints 
and connectors using microwave reflectometer techniques. 

R. W. Beattv, G. F. Engen, and W. J. Anson, IRE Trans 
Instrumentation 1-9, 219-226 (Sept. 1960). 
The reflection and loss of a waveguide joint are quantified 
by its VSWR and efficiency, respectively. These are con- 
veniently and accurately determined by techniques which 
employ a reflectometer with auxiliary tuners and an adjustable 
sliding termination. Depending upon the stability and gain 
of the associated apparatus, measurements of VSWR's of 1.001 
and lower, and of efficiencies of 99.99 per cent and higher 
may be obtained. The attenuations of short lengths of wave- 
guide may be determined at the same time that the efficiency 
measurement is made. 

A precision RF attenuation calibration system, C. M. Allred 
and C. C. Cook, IRE Trans. Instrumentation 1-9, 268-274 
(Sept. 1960). 

A new precision attenuation calibrating system with greatly 
increased sensitivity, stability, and measurement range has 
been completed at the National Bureau of Standards, Boulder 
Laboratories. The increased stability and sensitivity are 
achieved by the use of a highly accurate piston attenuator 
and precision phase shifter combined into a null system. 
The extended attenuation measurement range has been ob- 



tained by using a new mode launching system which is excited 
by essentially a constant current source of very high magni- 
tude. The padding necessary for correct impedance matching 
is kept at a minimum by use of a special noninteractive com- 
bining network. 

A precision RF power transfer standard, P. A. Hudson, 
IRE Trans. Instrumentation 1-9, 280-283 (Sept. 1960). 
A fast-reading, precision RF power transfer standard is 
described which is usable at fixed frequencies in the range 
10 Mc to 1000 Mc. The dynamic range is 1 watt to 1000 
watts with a resolution of 0.1 per cent to 0.01 per cent. 
Stability is 0.1 per cent per year. The standard consists of 
three directional couplers to cover the frequency range and 
a series of vacuum thermoelement detectors matched to 50 
ohms at each frequency interest. The thermoelements are 
connected to the secondary line of the appropriate coupler 
and provide a dc output voltage which is calibrated in terms 
of RF power in the coupler main line. 

A series of small, portable wattmeters employing vacuum 
thermoelements are also described. 
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R. S. Roth, J. Am. Ceram. Soc. 44, No. 1, 49-50 (Jan. 1 961). 
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Hydrogen-bonding in calcium-deficient hydroxyapatites, A. S. 

Posner, J. M. Stutman, and E. R. Lippincott, Nature 188, 

No. 4749, 486-487 (Nov. 1960). 
Hydrogen sulfide precipitation of the elements from 0.2-0.5 

normal hydrochloric acid, J. I. Hoffman, Chemist Analvst 
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Dielectric constant and dielectric loss of Ti0 2 (Rutile) at low 

frequencies, R. A. Parker and J. H. Wasilik, Phvs. Rev. 

120, No. 5, 1631-1637 (Dec, 1960). 
Generating functions for formal power series in non-commut- 
ing variables, K. Goldberg, Proc. Am. Math. Soc. 11, No. 

6, 988-991 (Dec. 1960). 
Tests for regression coefficients when errors are correlated, 

M. M. Siddiqui, Ann. Math. Stat. 31, No. 4, 929-938 

(Dec. 1960). 
Application of the method of polarized orbitals to the scat- 
tering of electrons from hydrogen, A. Temkin and J. C. 

Lamkin, Phys. Rev. 121, No. 3, 788-794 (Feb. 1961). 
A recommended standard resistor-noise test system, G. T. 

Conrad, Jr., N. Newman, and A. P. Stansbury, IRE Trans. 

Component Parts CP-7, No. 3, 71-88 (Sept. 1960). 
On the absorption spectrum of CF 2 and its vibrational 

analvsis, D. E. Mann and B. A. Thrush, J. Chem. Phys. 

33, No. 6, 1732-1734 (Dec. 1960). 
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